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Abstract 

We review the basic ideas and results on the vortex dynamics in clean 
superfluid Fermi systems. The forces acting on moving vortices are dis- 
cussed including the problem of the transverse force which was a matter of 
confusion for quite a time. Finally, we formulate the equations of the vor- 
tex dynamics which include all the forces and the inertial term associated 
with excitations bound to the moving vortex. 

1 Introduction 

Study of the vortex dynamics opens promising perspectives for understanding 
the most fundamental properties of condensed matter especially for supercon- 
ductors and superfluids. Among them, clean systems, i.e., the systems where 
the mean free path of quasiparticles is much longer than the characteristic co- 
herence length, offer more intriguing physics than dirty superconductors. For 
example, one of the fundamental problems can be formulated as follows: Speak- 
ing of clean superconductors one can, in particular, think of such a system where 
no relaxation processes are available, i.e., the mean free path of excitations is 
infinite. In particular, vortices in such systems should move without dissipation 
since there is no mechanism to absorb the energy. The vortex velocity should 
then be parallel to the transport current, which makes the induced electric field 
perpendicular to the current; the dissipation thus vanishes j • E = 0. This con- 
trasts with what we know of dirty superconductors: the vortex motion there is 
dissipative, and each moving vortex experiences a large friction; it generates an 
electric field parallel to the transport current which produces energy dissipation. 
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Clearly, a crossover should occur from dissipative to non-dissipative vortex mo- 
tion as the quasiparticle mean free path increases. The question is what is the 
condition which controls the crossover? 

This question is of a fundamental importance for our understanding of the 
dynamics of superconductors and, in a more broad sense, for the understanding 
of dynamic properties of quantum condensed matter in general. To illustrate 
the problem one can consider a simple example as follows. One can argue that 
a time dependent non-dissipative superconducting state, similarly to any other 
quantum state, can be described by a Hamiltonian dynamics based on a time 
dependent Schrodinger equation. Such a description for a weakly interacting 
Bose gas has been suggested by Pitaevskii and Gross Jl], |). It is widely used 
for superfluid helium II as well. The Gross-Pitaevskii equation is essentially a 
nonlinear Schrodinger equation, it has the imaginary factor — ih in front of the 
time derivative of the condensate wave function d^>/dt. On the other hand, the 
time dependent Ginzburg-Landau model which is a particular case of a more 
general Model F dynamics j| is believed to describe a relaxation dynamics of 
superconductors near the transition temperature. In contrast to the Gross- 
Pitaevskii equation, it has the time derivative /dt with a real factor in front 
of it. The question which we are interested in can be formulated as follows: 
What is the condition when the imaginary prefactor transforms into a real one? 

It seems that there is no universal answer to this simple question in general. 
However, the problem of crossover from non-dissipative to dissipative behavior 
of a condensed matter state can be solved for the particular example of vortex 
dynamics. It is known that the relaxation constant in the time dependent 
Ginzburg-Landau model has in fact a small imaginary part |5|, which 
results in appearance of a small transverse component of the electric field with 
respect to the current. We shall see later that the transverse component of 
the electric field increases at the expense of the longitudinal component as the 
mean free path of excitations grows. The crossover condition, however, does 
not coincide simply with the condition which divides superconductors between 
dirty and clean ones. The criterion rather involves the spectrum of excitations 
localized in the vortex cores; the distance between their levels takes the part of 
the energy gap. The condition for a non-dissipative vortex motion requires that 
the relaxation rate of localized excitations is smaller than the distance between 
the levels. This implies a much longer mean free path of excitations than the 
condition for a superconductor to be just in a clean limit. 

In the present paper we review the basic ideas and results on the vortex 
dynamics in clean Fermi systems. We concentrate largely on superconductors, 
the case of superfluid 3 He can be easily incorporated if one takes the limit of 
zero charge of carriers in the final results (to be discussed it in more detail 
later). We consider forces acting on moving vortices and clarify the conditions 
for the crossover from dissipative to non-dissipative vortex dynamics. A great 
deal of attention is given to the problem of the so called transverse force, i.e., 
the force perpendicular to the vortex velocity, which was a matter of confusion 
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for quite a time since its discovery. Finally, we formulate the equations of the 
vortex dynamics which include all the forces and the inertial term associated 
with excitations bound to the moving vortex. 

2 Boltzmann kinetic equation approach 

The forces on a vortex come from several different sources including the hydro- 
dynamic Magnus force, the force produced by excitations scattered from the 
vortex, and the force associated with the momentum flow from the heat bath 
to the vortex through the localized excitations. The whole rich and exciting 
physics involved in the vortex dynamics can be successfully described by the 
general microscopic time dependent theory. Here we rather discuss a general 
picture using a simple semi-classical approach based on the Boltzmann kinetic 
equation. The semi-classical approach assumes that the wavelength of exci- 
tations is much shorter than the superfluid coherence length, pp£ 3> 1 (from 
now on we put % = 1). This condition is quite safely fulfilled in almost all 
superconductors and in superfluid 3 He. Only in some high-temperature super- 
conductors, its accuracy may be not very good. For simplicity, we consider 
s-wave superconductors. We concentrate on isolated vortices such that their 
cores do not overlap, i.e., on the region of magnetic fields H <C H C 2- 

One can distinguish two types of excitations: Excitations localized in the 
vortex core, and those which are not localized but move in the vortex potential 
under the action of magnetic filed. We start our discussion with the localized 
excitations. 

2.1 Localized excitations 

We remind that the profile of the order parameter A (r) near the vortex core 
produces a potential well where localized states with a discrete spectrum exist. 
The localized states correspond to energies |e| < Aoo. The spectrum has the so 
called anomalous branch with the radial quantum number n — whose energy 
varies from -A M to +Aoo as the particle impact parameter b changes from — oo 
to +oo and crosses e = 0, being an odd function of b. For low e <C A^, the 
anomalous branch is £q — —loo^ where /i = —bp± is the angular momentum, and 
p± is the momentum in the plane perpendicular to the vortex axis (see also 
Fig. [l]). In a s-wave superconductor with an axisymmetric vortex, the angular 
momentum \x is quantized and so is the spectrum, loq being the distance between 
the discrete levels in the vortex core. The spectrum also has branches with n =/= 
which are separated from the one with n — by energies of the order of A. They 
are even in b and symmetric in energy with respect to e = 0. We denote the 
separation between the levels with neighboring angular momenta through 



A 

8 




Figure 1: Spectrum of excitations localized in the vortex core. 



The interlevel spacing co n (b) is an even function of b for n = 0, and it is an odd 
function for n =/= 0. 

We choose the direction of the z axis in such a way that the vortex has a 
positive circulation. The z axis is thus parallel to the magnetic field for positive 
charge of carriers, and it is antiparallel to it for negative charge: z = hsign (e). 
Since the particle velocity vj_ in the plane perpendicular to the vortex axis 
makes an angle a with the x axis, the cylindrical coordinates of the position 
point (p, <j>) are connected with the impact parameter and the coordinate along 
the trajectory through p 2 = b 2 + s 2 where 

b = psm(<fi — a) ; s = pcos((f> — a). (2) 

The coordinates are shown in Fig. ||. 

The first step is as follows. We assume that the quasiclassical spectrum e n (b) 
of a particle plays the role of its effective Hamiltonian. We can thus invoke the 
Boltzmann equation in the canonical form 



m df de n de n df ( df 
dt da dp da dp \ dt 



(3) 

coll 



to describe the quasiparticle distribution B. Equation (0) has been derived in 
Ref. |9) from the set of microscopic kinetic equations. 

In the time derivative, the energy e„ contains a time dependence through 
p{t) = [(r — VLt) x p] • z such that 

0f(0) df^de n dp dfWde n . { 

~dT = ~dT~d~iI~dt = ~dT~diI ([Pi x Vi] ' z) • (4) 
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Figure 2: Transformation from cylindrical coordinates to the frame associated 
with a particle moving in the vortex core. The line AB is the particle trajectory 
making an angle a with the x axis and passing through the position point (p, <ft) 
at a distance b from the vortex axis. 



Kinetic equation takes the form 

dft 9f(°) de n dh de n de n df t fdf\ 

where the distribution function is separated into an equilibrium and a non- 
equilibrium parts, f = f + f*i, respectively. Here = 1 — 2n e = tanh (e/2T) 
with n e being the Fermi function. 

We shall simplify the collision integral in Eq. (^) using the relaxation-time 
approximation 

oil Tn 



dt 



where r„ ~ r. With this approximation, the mean free time can be of any origin. 
For definiteness, we assume that the most effective relaxation is brought about 
by impurities, as is the case in almost all practical superconducting compounds. 
However, the relaxation time r can be due to the electron-phonon interaction, 
as well, or due to quasiparticle-quasiparticle scattering. The latter is, in fact, 



the only relaxation mechanism available in 3 He |10| . Therefore, the term "clean" 
does not necessarily mean a low concentration of impurities but simply refers 
to a situation when the mean free time is long. 

Equation (||) with the collision integral in the form of Eq. (||) is easy to 
solve. For an axi-symmctric s-wave vortex the energies e„ do not depend on 
a and the term de n /da vanishes. Let us take the distribution function in the 
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form 

h = g— [70 ([vi x pj-z) +7h(v l • p_l)J (7) 

where the factors 70, h are to be found from the Boltzmann equation (^|). The 
result for a steady vortex motion is 

T° = 2 "2 7 1 ' ^ H = 2 ? 2 " 1 ■ ( 8 ) 

This generalizes the result first obtained in Rcf. jllj]. 
2.2 Delocalized excitations 

A delocalized particle moves mostly far from the vortex core where the order 
parameter is constant and the superfluid velocity potential pv s is small com- 
pared to A. Kinetic equation for delocalized excitations can thus be written as 
for a particle in a magnetic field with a semi classical spectrum 



e P + (9) 

where £ p = p 2 /2m — Ep. As shown in Ref. iQ the kinetic equation for a 
particle moving in a vortex array has the conventional form 

dt_ df_ v df_ _ (M_ 
dt dp 9 dr \dt 

The force is the elementary Lorentz force 

dp e uj c 
f = — = -v g x H =— p F x z (10) 
at c g 



where oj c — |e| H/mc is the cyclotron frequency, 

de p _ v F 



dp g 

is the group velocity, and 

e 

The driving term can be written as 

df _ dfW de 
dt ~ de dt 

where 



(11) 



de e rTT . lu c 

— = ev g ■ E =-v g ■ H x v L = — p ± ■ z x v L 

dt c g 



G 



The kinetic equation becomes 

df,df^Lo c dfi uj fdf\ 

■nr + -5 Pi • z x vi, + ^ \p F xz]=[ — ) . 12 

We omit the spatial derivative of the distribution function since it is constant 
in space. 

For the energy spectrum of Eq. (61) the collision integral is O] 

I) =-"*■ 

, 5 Vcoii 9r 
Kinetic equation (|l2] ) takes the final form 



.ef<°> 54 1 . 

p_L • Z X V L \ — — = fi. 

oe da uj c t 



Its solution is Eq. (Q) with 



TO = 2 2,1 ' = 22,1 - ( 13 ) 



3 Forces 



At the second step we calculate the force acting on a vortex from the environ- 
ment. This force is exerted via excitations which travel near the vortex through 
transfer of their momenta to the vortex. Consider first localized excitations. 
The transferred momentum is 

F (ioc) = 1 >P f d P z dad ^ dP" , f dpz dadfl d€n r - If (~\A\ 

2 2^ J 2tt 2tt dt 1 2 <W 2tt 2^ 96 [Z V±J X ' 1 j 

n n 

Here we make use of the Hamilton equation 
dPn _ de n 

¥ = - Ve " = -¥ |zxVi| - (15) 

The second equality follows from the fact that, in the coordinate frame (s, b) of 
Fig. the energy only depends on the particle impact parameter b. 

The transfer of the momentum from delocalized excitations in the semi- 
classical approximation has the form 

F (dci) = I V f dPz ( — ^ t 



2 ^ J 2tt / 2tt dt 

n 

where dp n dt is given by Eq. (jTc|). The force becomes 
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The sum is over quantized states of a semi-classical quasiparticle in a magnetic 
field. Indeed, the energy spectrum Eq. (^|) becomes quantized (l2) in a magnetic 
field and turns into a modified Landau spectrum e n such that 



= io c n - p]_/2m. 



The quantity 



9_ 



On 
d7 



is thus the density of states. As a result, 

E- 

n 

and the force becomes [Ol 



dn , 
Te de 



F (dol) 
cnv 

The total force is 



e>A 



dp z 
2tt 



/g[ Pj _xz]f, 



p(loc) i p(dcl) 



(16) 
(17) 



With the Ansatz (^), the total force Eq. ( |l7j ) splits into two terms F e 



with the friction Fn and transverse Fj_ forces given by 



-ttN 



E 



df (°) d M \\ A . A 
Wn70 — —)>+(!- tanh 



^o7h 



<9e 2 



2T 



y 



vl,(18) 



A 



9e Ty/M 1 "^^ J7h 



where N is the quasiparticle (electron) density, 
Fermi surface with the weight 7rp^, 



= Vp 1 J np\d Pz (...), 



[zxv L ], (19) 
is the average over the 

(20) 



and Vf is the volume encompassed by the Fermi surface. For an isotropic Fermi 
surface, 



((...)) = I J S m*9dd (...). 



where p z = pf cos 6. Only the spectral branch with n = contributes to the 
transverse force Fj_ as all w„ with n ^ are odd functions of \x and thus drop 
out of the sum over n in Eq. (^) . 

We emphasize that the force F on v is defined as the response of the whole 
environment to the vortex displacement. It is therefore the total force acting on 
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the vortex from the ambient system, including all partial forces such as the lon- 
gitudinal friction force and the non-dissipative transverse force. The transverse 
force, in turn, includes various parts which can be identified historically [ fL5| as 
the Iordanskii force the spectral flow force [|l7|, [l^, |l9|, and the Magnus 
force. We shall discuss this later in more detail in SectioirT 



3.1 Flux flow conductivity 

The force from the environment is balanced by the Lorentz force: 

F L = -^[j tr xz]sign(e), (21) 

with the flux quantum $o = 7re/|e|. The force balance equation 

F L + F cnv = (22) 

determines the transport current in terms of the vortex velocity and thus allows 
us to find the flux flow conductivity tensor. The longitudinal force F|| defines the 
friction coefficient in the vortex equation of motion and determines the Ohmic 
component of the conductivity o~o- Expressing the vortex velocity Vi through 
the average electric field E, as Vi = c [E x i}/B sign(e), we find 



N\e\c 
B 




1 — tanh 



Aqq 

2T 



LO c T 



The transverse force determines the Hall conductivity 



oh 



Nec 



,2^2 



1 de 



de 
~2~ 



1 — tanh ■ 



(23) 



(24) 



The main conclusion is that the Ohmic and Hall conductivities depend on 
the purity of the sample through the parameters ujqt (remind that uj n ~ ujq) 
and io c t. Note that lo c ~ (H j ' H C 2)ujq thus lj c <C t^o- One can distinguish two 
regimes: moderately clean ujqt 1 and superclean loqt 1. Note that the 
moderately clean regime still requires that the superconductor is clean in the 
usual sense A^t 3> 1. 

In the moderately clean limit where ujqt 1, the conductivity roughly 
follows the Bardeen and Stephen expression [^0) at low temperatures though it 
exhibits an extra temperature-dependent factor Aqo/Tc on approaching T c [ pTj 



H C 2 



The factor A^/Tc appears because the number of delocalized quasiparticlcs 
contributing to the vortex dynamics decreases near T c . This extra factor has 
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been recently identified experimentally [Ell. The Hall conductivity and the 
Hall angle are small cth ~ (ojot) <jo and tanO^r ~ (wot), respectively. The 
contribution from delocalized states are not important since lu c < wo- 

In the superclean limit luqt 3> 1, on the contrary, the Ohmic conductivity is 
small. The vortex dynamics becomes non-dissipative. In particular, if uj c t <C 1, 
the corresponding Hall conductivity is 

Nee. 

an = tanh — — . 

B \2T J 

If lo c t 3> 1, the hyperbolic tangent should be replaced with unity. 



4 Transverse force 

Let us now discuss the forces acting on a moving vortex in more detail. The 
friction force is determined by the Eq. ([l8]) . It is proportional to the mean free 
path of excitations in the moderately clean regime and vanishes in the superclean 
limit. The transverse force Eq. (|l9| ) deserves a more careful discussion because 
it has been a matter of controversy for a long time since first calculated for 
vortices in helium II by Lifshitz and Pitaevskii |22) and then by Iordanskii Jl6| . 
The review |Q tells about old disputes (see also |23|). Recently, the presence 
of the transverse force has been questioned in Refs. 



mm- 

As we see, however, the microscopic picture gives a finite transverse force in 
a full accordance with the symmetry arguments which allow a transverse force 
in a chiral system such as a moving vortex in presence of a superflow. Of course, 
the magnitude of the transverse force is not an universal quantity. It appears 
to depend on parameters of the system. Indeed, in the superclean limit, when 
both lu c t and ojqt are much larger than unity, the factors 7h — 7h — 1; and the 
transverse force Eq. (|l^) becomes 

F_l = nN[i x Vi ]. (25) 

The balance Eq. (j^) gives the transport current in the form 

j tr = Nev L . (26) 

This equation is consistent with the Hclmholtz theorem of conservation of cir- 
culation in an ideal fluid: vortices move together with the flow. The correction 
to the distribution function of excitations is simply 

fi = • px) (27) 

such that the full distribution function is f = (e — • p): excitations move 
together with the vortex being in equilibrium with the vortex array. To under- 
stand this let us consider first the mean free path of delocalized excitations with 
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respect to their collisions with vortices. If the vortex cross section is a v , the 
mean free path is £ v — l/a v n v where n v — i?/$o is the density of vortices. We 
shall see in a moment that the vortex cross section is a v ~ p F x so that 

tv ~ Pf/ti v ~ v f /lu C i (28) 

i.e., t ~ tl where tl is the Larmor radius. In the limit uj c t 3> 1, the vortex 
mean free path £ v becomes shorter than the impurity mean free path £i mp = vpT 
so that the delocalized excitations scatter on vortices more frequently than on 
impurities and thus come to equilibrium with moving vortices. A similar con- 
sideration also applies to localized excitations: In the limit loqt 3> 1 interaction 
with a vortex is more effective than relaxation on impurities, and the excitations 
come to equilibrium with the moving vortex. 

In the intermediate case when oj c t — > but luqt ~ 1, only the localized 
excitations are out of equilibrium. The excitations with energies above the 
energy gap Aoo have £ v £- lmp and are now in equilibrium with the heat bath: 
both 7q = 0, and 7^ =0. As a result these excitations do not influence 
considerably the vortex motion. Delocalized excitations, however, do affect the 
vortex motion if uj c t is comparable with unity. Finally, in the limit u>qt, uj c t <C 1 
the transverse force Eq. ( |l9] ) disappears. The longitudinal force dominates: the 
vortex dynamics is dissipative. 

In general, we can identify several contribution to the transverse force. One 
can present the full transverse force Eq. (fL9|) in the form 



The force 



-nN a [v L xz]-7rA^„[v L xz] + F sf . (29) 



+-kn(\- tanh ^) -pL^ [v L x z] (30) 

is called the spectral flow force Jl7| . 

The spectral flow force F s f is due to the momentum flow from the Fermi sea 
of normal excitations to the moving vortex via the gapless spectral branch going 
through the vortex core from negative to positive energies |l8),[^9|. Due to the 
time dependent angular momentum of the excitations /j,(t) = [(r — v^<) x p] • z, 
there appears a flow (with the velocity dfi/dt) of spectral levels characterized 
by the angular momentum /i. Each particle on a level carries a momentum pp. 
The momentum transfer to the vortex is effective if the quasiparticle relaxation 
occurs quickly: the factor {uJqTq + 1) _1 in the first line of Eq. ( |30| ) accounts 
for the relaxation on localized levels: the relaxation and hence the momentum 
transfer is complete for luqt <C 1, and vanishes in the opposite limit. The first 
term in Eq. (p0|) thus describes a disorder-mediated momentum flow along 
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the anomalous chiral branch Eq (fi) for energies below the gap. The second 
line in Eq. (^) accounts for the spectral flow for energies above the gap. The 
corresponding factor which takes into account the relaxation rate is (u>;?t 2 + 1) _1 . 

Note that F onv is calculated in the reference frame where the normal com- 
ponent is at rest v„ = for stationary vortices. In superconductors, the normal 
velocity is always zero in equilibrium, the excitations being at rest in the refer- 
ence frame associated with the crystal lattice. In superfluid 3 He, the situation is 
similar: the normal component has a rather large viscosity so that it is normally 
at rest in the container frame of reference. For a nonzero v„, we would need to 
replace V/, with vj, — v„ everywhere in F cnv . Simultaneously, the Lorentz force 
would also include the contribution due to the normal velocity. 

Let us now turn to the force balance equation (^|]) . In presence of an electric 
field, the transport current is not entirely due to a supercurrent, a part of it being 
carried by delocalized quasiparticles. Far from the vortex core, the quasiparticle 
current is 

j(qp) = _ 2v (0) e / V ± . 9 fid£^ (31) 

where ^(0) is the single-spin normal density of states, and d£l p is the elementary 
solid angle in the direction of the momentum p. Using Eqs. (0) and we 
find 

jfori = N n e^v L + N n ei Q [z x y L ] (32) 
where the density of normal quasiparticles is 

N n = N I g— de. (33) 



e>A 



Writing the transport current as jt r = N s ev s +j( qp ) we get the force balance 
Eq. (p2|) in the form 



Here 

is the Magnus force; 



^ M + F[ qp) + F I + F 8f + F,| = 0. (34) 



F M = ttN s [(v g - v L ) x z] 



F (qp) = £0 [j(q P ) x £] sign (e) (3 5 ) 

is the Lorentz force from the quasiparticle current Eq. (|3^), and 

Fi = nN n [(v n - vx) x z] 

is called the Iordanskii force [^6| . The Iordanskii force is the counterpart of the 
Magnus force for normal excitations. We have included the normal velocity v n 
to make this similarity more transparent. 
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The spectral flow force vanishes in the limit r — > oo when both loot 3> 1 
and lu c t ^> 1. In this limit the transverse force is given by Eq. (p5f). The 
quasiparticle current is j' qp * ) = N n evL so that the force from the quasiparticle 
current compensates the Iordanskii force. The force balance Eq. (|?4|) reduces 
to Fm = 0. The vortex thus moves with the superfluid velocity vj, = v s . As 
follows from Eq. (^?]), quasiparticles also have a velocity so that all the 
particles move together which amounts to the total current as in Eq. (|26|). 

On the contrary, the spectral flow force has its maximum value for moder- 
ately clean limit, loqt -C 1. Equation ( |30| ) gives in this limit 

F sf = nN[v L x z]. 

This completely compensates the first two terms in Eq. (p9|), i.e, the Iordanskii 
force and the part of the Magnus force that contains the vortex velocity; the 
transverse force vanishes. The quasiparticle current vanishes even faster because 
lu c <C u>o- The Lorentz force is balanced only by a friction force F||. As a result, 
the dissipative dynamics is restored. 



4.1 Low-field limit and superfluid 3 He 

The low-field limit when lu c t <C 1 is most practical for superconductors. More- 
over, this regime is realized in electrically neutral superfluids such as 3 He. At 
the first glance, it is simply because u> c vanishes together with the charge of 
carriers. However, this is not completely correct. In fact, to estimate a de- 
viation from equilibrium of dclocalizcd excitations in this case one has again 
to compare the mean free path of excitations with their mean free path with 
respect to scattering by vortices. Keeping in mind that the vortex density is 
n v = 20,/ k where Q is an angular velocity of a rotating container and k = n/m 
is the circulation quantum, Eq. ( p8| ) gives £ v ~ vp/Vt. We observe that the 
cyclotron frequency is replaced with the rotation velocity in a full compliance 
with the Larmor theorem. The ratio of the particle-particle mean free path £ 
to the vortex mean free path is £/l v ~ Or. With the practical rotation veloc- 
ity £1 of a few radians per second one always has £ v exceedingly larger than I, 
Delocalized excitations are thus at rest in the container frame. 

Consider this regime in more detail. Since delocalized excitations are in 
equilibrium the quasiparticle current vanishes. The force balance becomes 

F M + Fi + F sf + F,| = (36) 

where the spectral flow force is 

+ttN M - tanh ^ J [v L x z] . (37) 
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It is interesting to observe that the spectral flow force from delocalized states 
in this case is related to the anomalous contribution to the transverse vortex 
cross section for scattering of delocalized quasiparticles. The vortex cross sec- 
tions were calculated in Refs. [(26l Wf\. The transverse cross section is 



7T 

(T_L = 

P± 



(38) 



per unit vortex length. Note that the transport cross section which is responsible 
for the scattering contribution to the longitudinal force is vanishingly small in 
the semi-classical limit. Inserting Eq. ([38|) into the expression for the force 
exerted on the vortex by scattered excitations, 

, (sc) _ f <kdf<® fd Pz 3 



we recover those contributions to Eq. ( p9| ) which are due to the normal exci- 
tations, namely the term —irN n [v^ x z] and the last term in the spectral flow 
force Eq. @. 

The first term in Eq. (^8|) corresponds to the cross section of a vortex in a 
Bose superfluid (l5) 

<7j_ = 2-K/m B v g (40) 

where v g is the group velocity which, in our case, is defined by Eq. (|ll|), and 
tub is the mass of a Bosonic atom. Note that, in our case, tub = 2m. The 
corresponding part of the transverse cross section can be easily obtained from 
the semi-classical description. Indeed, the Doppler energy due to the vortex 
velocity is p • v s = p • \7x/ m B- Its contribution to the quasiparticle action is 

A = -±[ v .V x dt = —^fgd, = —>-8 X . 
tub J m B v g J ds m B v g 

Here s is the coordinate along the particle trajectory as in Fig. |[ and S\ is 
the variation of the order parameter phase along the trajectory. The change in 
the transverse momentum of the particle is Sp± = dA/db hence the transverse 
cross section becomes 



>± 



6 J±db= A+ - A - 
P P 



where A± is the action along the trajectory passing on the left (right) side of 
the vortex. Since 5x+ — = ~27r we recover Eq. ([io|). With this expression 
for the cross section, Eq. ( p9|) g ives the Iordanskii force Fi = — TrN n [vL x z]. 

The second term in Eq. (|38D originates |2(| from the fact that here, as dis- 
tinct from the situation in a Bose superfluid, the phase of the single-particle 
wave function changes by it upon encircling the vortex, while it is the order 
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parameter phase which changes by 2tt. It is this singularity produced by the 
vortex in the single-particle wave function which results in the anomalous con- 
tribution to the cross section in Eq. (|38|). Inserted into Eq. (|39|), it exactly 
reproduces the second term in Eq. (|37j). We see that the spectral flow force is 
related to a single-particle anomaly associated with the vortex. 



5 Vortex momentum 

The vortex mass in superfluids and superconductors has been a long standing 
problem in vortex physics and remains to be an issue of controversies. There are 
different approaches to its definition. In early works on this subject, the vortex 
mass was determined through an increase in the free energy of a superconductor 
calculated as an expansion in slow time derivatives of the order parameter. The 
quasiparticle distribution was assumed to be essentially as in equilibrium. First 
used by Suhl |^8| (see also [^9| ) this approach yields the mass of the order of one 
quasiparticle mass (electron, in case of superconductor) per atomic layer. An- 
other approach consists in calculating an electromagnetic energy E 2 /8ir which 
is proportional to the square of the vortex velocity. This gives rise to the so 
called electromagnetic mass which, in good metals, is of the same order of 
magnitude (see Ref. |5l) for a review). 

A crucial disadvantage of the above definitions of the vortex mass is that they 
do not take into account the kinetics of excitations disturbed by a moving vortex. 
We shall see that the inertia of excitations contributes much more to the vortex 
mass than what the old calculations predict. The kinetic equation approach 
described here is able to incorporate this effect. To implement this method 
we find the force necessary to support an unsteady vortex motion. Identifying 
then the contribution to the force proportional to the vortex acceleration, one 
defines the vortex mass as a coefficient of proportionality. This method was 
first applied for vortices in superclean superconductors in Ref. |32| and then 



was used by other authors (see for example, Refs. [|33, 34 ). The resulting mass 
is of the order of the total mass of all electrons within the area occupied by 
the vortex core. We will refer to this mass as to the dynamic mass. Since 
the dynamic mass originates from the inertia of excitations localized in the 
vortex core it can also be calculated through the momentum carried by localized 
excitations |?5| . We shall see that dynamic mass displays a nontrivial feature: 
it is a tensor whose components depend on the quasiparticle mean free time. 
In s-wave superconductors, this tensor is diagonal in the superclean limit. The 
diagonal mass decreases rapidly as a function of the mean free time, and the off- 
diagonal components dominate in the moderately clean regime. These results 
were obtained in Ref. |3(| . Our results agree with the previous work |33|, Q 
in the limit r — > oo. 
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5.1 Equation of vortex dynamics 

To introduce the vortex momentum we consider a non-steady motion of a vortex 
such that its acceleration is small. We again start with the delocalized excita- 
tions. Multiplying Eq. (|^) by p_i_/2 and summing up over all the quantum 
numbers, we obtain 

F£ v c) = Fff - d -^ r nN tanh ^ j [v L x z] (41) 

where the l.h.s. of Eq. (^) is the force from the environment on a moving 
vortex Eq. @. 

The first term in the r.h.s. of Eq. j4l| ) is the force exerted on the vortex by 
the heat bath via excitations localized in the vortex core: 

p(ioc)__l\^ f n ( df\ dp z dadji 



■ coll 



1 X^f n f df \ d Pz dadv 



The second term in the r.h.s. of ( |4l] ) is the change in the vortex momentum 

r<^-I E / P±/l «. ,43) 

n 

We turn now to delocalized excitations. Multiplying Eq. ( fl2| ) by pj_/2 again 
and summing over all the states we find 

F (dd) F (dei) 9P( del ) l v [dp z dadfWu> 

n nv j - Peon -—qt + ^J -^2^^j p±{p± - [zxvl]) 



1 — tanh 



Aoo 

2T 



[v L xz]. (44) 



Here 



-(del) 1 /" dp z da /" fi dp z da 

u = -2 ? y p - UL ^ ^ = / >A _ ^ ^ 

is the force from the heat bath, and the corresponding contribution to the vortex 
momentum is 



,(dcl) 



= --T [p±h^- = -[ p^h^-de. (45) 
2^ J V±J1 2tt 2tt Je>A oa w c 2tt 2tt V ' 



The total momentum is P = P< loc ) + p( dcl ). 

After a little of algebra, the total force from the heat bath F co n = F^°^ + 



PcoU can ^ e written in the form 



Faf + F|| 
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where the friction and spectral flow forces are determined by Eqs. (jig) and 



30|), respectively. The total force from environment Eqs. (41) and ( [44[ ) takes 
the form 

F nv = F rf + Fy - — - tt7V[v l x z] . (46) 

This equation agrees with Eq. ( p9| ) for a steady motion of vortices. 

The equation of vortex dynamics is obtained by variation of the supercon- 
ducting free energy plus the external held energy with respect to the vortex 
displacement. The variation of the superfluid free energy gives the force from 
the environment, F onv , while the variation of the external field energy produces 
the external Lorentz force. In the absence of pinning the total energy is trans- 
lationally invariant. Therefore, the requirement of zero variation of the free 
energy gives again the condition Fl + F env = in the form the force balance. 
Using our expression for F C11V , the force balance can now be written in the form 
similar to Eq. (p|) 



•M+F^ p) +F I + F sf + F|| =^ (-17) 



where the r.h.s. contains the time derivative of the vortex momentum due to 
a non-steady vortex motion. For a steady vortex motion, Eq. (WTn reduces to 



Eq. (p4[). The physical meaning of the Eq. (47) is simple. The l.h.s. of this 
equation accounts for all the forces acting on a moving straight vortex line. The 
r.h.s. of Eq. ( fi7j ) comes from the inertia of excitations and is identified as the 
change in the vortex momentum. The definition of Eq. ( p3| ) is similar to that 
used in Refs. M, [3pJ. Note that the delocalized quasiparticles do not contribute 
to the vortex momentum in the limit lu c t <C 1 because they are in equilibrium 
with the heat bath. 



5.2 Vortex mass 

Having defined the vortex momentum, we calculate the vortex mass. The 
distribution function is given by Eq. (0). The vortex momentum becomes 
Pi = MikUk] it has both longitudinal and transverse components with respect 
to the vortex velocity. For a vortex with the symmetry not less than the four- 
fold, the effective mass tensor per unit length is = M\\8ik — M±eikjZj where 
M|| = M|| e + Mn h and M± = M± e — M±f l - Each component contains contri- 
butions from localized and delocalized states such that M\\ t ± — M9°^ + M, 1 ^ 1 ' 
where 

KS — ((? / ™ ^ f » 

(48) 
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and 

<H/ Pli^^^d, (49) 

The same expression holds for M± Bi h where 7h is replaced with 70. The indexes 
e, h indicate the corresponding momentum integrations over the electron and 
hole parts of the Fermi surface, respectively. Only the branch with n — gives 
the contribution to the transverse mass because 70 is odd in fi for n^O. 

If the vortex acceleration is slow, one can use expressions Eqs. (||), ( fL3| ) 
for a steadily moving vortex to calculate the vortex inertia. Consider first the 
contribution of the states with |e| > A^. Since 7^ do not depend on energy 
and momentum, Eq. (^) gives 

M l ;r = ^^7H.o=^' i) «,o 

The contribution of the delocalized states decreases as lu c t gets smaller. In the 
limit of vanishing lo c t the vortex mass is determined by localized excitations. 
The localized excitations dominate also at low temperatures T <C Aoo . One has 
in this case 

M ]le , h = *N**> , M± e< h = nN^ 

For u) c t 3> 1 the mass tensor for delocalized states is diagonal Mjfc = MuSik 
where M|| dc ^ = mN n So', it is equal to the mass of normal particles in the area 
occupied by the vortex. The mass tensor for localized states becomes diagonal in 
the superclean limit where T 2 t/E f > 1 with M,| loc) - ttN ((c^ 1 )) - n^rnN. 
This is the mass of all electrons in the area occupied by the vortex core. The 
mass decreases with r. In the moderately clean regime T^t/Ef <S 1 where 
uj n T <S 1) the diagonal component vanishes as r 2 , and the mass tensor is domi- 
nated by the off-diagonal part. 

We should emphasize an important point that, in contrast to a conventional 
physical body, the mass of a vortex is not a constant quantity for a given system: 
it may depend on the frequency lo of the external drive. Indeed, for a nonzero 
u we find from Eq. (||) 

CJ 2 T 2 LO n T n {\ ~ lUJT n ) 

7H = n n , : To! 70 _ 



(1 - iUT n ) 2 ' uj^tI + (1 - iLOT n ) 2 ' 

As a result, all the dynamic characteristics of vortices including the conductivity 
and the effective mass acquire a frequency dispersion. In the limit t — » 00, poles 
in 70, h appear at a frequency equal to the energy spacing between the quasi- 
particle states in the vortex core which gives rise to resonances in absorption of 
an external electromagnetic field p2, p7[. 
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6 Conclusions 



All the dynamic characteristics of vortices such as vortex friction, the flux flow 
conductivity, the Hall effect, and the vortex mass are determined by the mean 
free path of excitations which interact with vortices. The key parameter is uqt 
where loq is the interlevel spacing for the quasiparticle states in the vortex core. 
For small uj t, vortices experience viscous flow; the Lorentz force is opposed 
by a friction force while the transverse force vanishes. On the contrary, in a 
superclean regime when uoqt » 1, vortices move with a superfiow as in an ideal 
fluid; the Hall angle is 7r/2, the friction force is zero, while the transverse force 
reaches its maximum value. The vortex mass is a tensor. The longitudinal 
component dominates in the superclean regime: it is the mass of all excitations 
in the vortex core. On the contrary, the transverse component is the largest one 
in the limit lo$t <C 1. 

I am pleased to acknowledge helpful discussions with E. Sonin and G. Volovik. 
This work was supported by the Russian Foundation for Basic Research (Grant 
No. 99-02-16043). 
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